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ABSTRACT 

We propose an alternative explanation for the post-breakout emission of SN 2008D associated with the 
X-ray transient 080109. Observations of this object show a very small contrast of 0.35 dex between 
the light-curve minimum occurring soon after the breakout, and the main luminosity peak that is 
due to radioactive heating of the ejecta. Hydrodynamical models show that the cooling of a shocked 
Wolf-Rayet star leads to a much greater difference (> 0.9 dex). Our proposed scenario is that of a 
jet produced during the explosion which deposits ^^Ni-rich material in the outer layers of the ejecta. 
The presence of high-velocity radioactive material allows us to reproduce the complete luminosity 
evolution of the object. Without outer ^^Ni it could be possible to reproduce the early emission 
purely from cooling of the shocked envelope by assuming a larger progenitor than a Wolf-Rayet star, 
but that would require an initial density structure significantly different from what is predicted by 
stellar evolution models. Analytic models of the cooling phase have been proposed reproduce the early 
emission of SN 2008D with an extended progenitor. However, we found that the models are valid only 
until 1.5 days after the explosion where only two data of SN 2008D are available. We also discuss 
the possibility of the interaction of the ejecta with a binary companion, based on published analytic 
expressions. However, the binary separation required to fit the early emission should be < 3i?0, 
which is too small for a system containing two massive stars. 

Subject headings: hydrodynamics — supernovae: general — supernovae: individual: SN 2008D 
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1. INTRODUCTION 

Supernova (SN) 2008D attracted a good deal of atten- 
tion because of its unusual observational characteristics. 
Most outstandingly, the serendipitous detection of the X- 
ray transient (XRT) 080109 associated with the SN ex- 
plosion, during the Swift follow-up o f SN 2007uy, another 
SN in the same galaxy, NGC 27 70 (iBerger fc Soderberd 
120081: iKong fc Maccarond [20081 ) . The optical counter- 
part of the XRT was revealed at the same position of the 
transient a few hours after, allowing to have unprece- 
dented early coverage of the SN emission. Initially the 
SN was classified as a broad- line Type Ic SN (SN Ic) 
(|Blondin et all l2008t iValenti et al] I2008D based on the 
broad absorptions and lack of hydrogen and helium in the 
first spectrum. Later on, however, the spectra revealed 
the presence of strong He I li nes, which changed the clas- 
sifica tion to Type lb (Modiaz et al.l[2008l : iMalesani et al.l 
MM)- 

Type lb SNe, as well as Type Ic and the transitional 
Type lib SNe, are believed to be the result of the core 
collapse of massive stars (A/zams ^ 10 Mq) that have 
lost most or all of their hydrogen (and often helium) 
layers before the explosion. F or that reason, they ar e 
called "stripped-envelope SNe" (|Clocchiatti et al.lll996f ). 
However, the mechanism by which the envelope is re- 
moved is not fully understood. Strong winds of sin- 
gle massive progenitors (M > 30 Mq), sudden erup- 
tions and binary interaction have been proposed as possi- 
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ble ex planation s (Macde r fc MevneSl 19941 : iNomoto et al.l 
[19951 IcTTOckTet al, ,20o¥ " among others), with the bi- 
nary origin being the most probable scenario. In all 
cases a Wolf-Rayet (WR) structure is expected before 
the explosion, although for SNe lib a thin H enve- 
lope is expect ed to remain before the explosion (e.g., 
see S N 1993J (INomoto et al.l 119931: [Podsiadlowski et alJ 
[1993), and SN 2011dh dBersten et al.l l2012t)V Interest 
in these objects has recently grown due to its connec- 
tion with long gamma-ray bursts (GRB). To date six 
SNe have been associated with GRB, one of which also 
showed XRT. In all of these cases a highly energetic Type 
Ic SN, also called "hypernova" was observed. Note, how- 
ever, that the term hypernova is more generally used, and 
refers to SNe with high luminosity or unusually broad 
lines, independently of their association with GRB. The 
origin of hypernovae is thought to be a rapidly rotatin g 
and accreting compact object |MacFadven et al.ll2001| ). 
or a magnetar ([Thompson et al.l2004l l. which should pro- 
duce relativistic outflows. 

The nature of the XRT associated with SN 2008D 
is controversial. Some authors are in favor of the su- 
pernova shock- breakout origin ([Soderberg et al.l [20081: 
iChevalier fc Fran sson 2008), while others consider that 
the transient was caused by a mildly relativistic jet pen- 
etrating through th e ei ivelope of the progen itor star 
(|Mazzah et al.l 120081: [lH [200a IXu et all I2008D . In the 
latter case, the XRT could represent a transition be- 
tween the most energetic hypernovae and standard core- 
collapse SNe. However, the presence of a GRB or 
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the radio detection of superluminal motions caused by 
a long-lived rel ativistic outflow were firmly ruled out 
for S N 2008D (jSoderberg et all [2008t iBietenholz eTall 
120091 ). Nevertheless, the lack o f confirmation of a 
thermal component for the XRT (jModjaz et all 120091 ) . 
and the s trong; evidence of a n asyrametric explosion 
(|Modiaz et al. 2009; Gorosabc fet all 120101), possibly 
bipolar ( Tanaka et al.i ■2009t .Maund et al.l 120091) . leaves 
room for the mildly relativistic-jet scenario. 

Another interesting feature of of SN 2008D was its 
double-peaked optical light curves (LC). The first peak 
occurred at sa 1 day after the XRT, and was similar to 
that of SN 2006a.i, w hich was connected with a GRB 
(jCampana et al.ll2006[ ). The second and main peak hap- 
pened at « 20 days after the XRT, consistently with 
other SNe lb and Ic, and its origin is related with the 
decay of radioactive material synthesized during the ex- 
plosion. By modeling the emission around the main peak, 
explosion parameters of SN 2008D, such as ejecta mass 
(Mej), kinetic energy (Ek), and nickel mass (M(^^Ni)) 
wer e estimated in previous studies. Using analytic mod- 
els, iSoderberg et all (|2008| ) found Moj = 3 - 5 M©, and 
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Monte Carlo simulations. iMazzali et a l. (2008) suggested 
that = 7 M©, Ek =6 foe, and Mi:'>^Wi) = 0.09 Mq. 
I Tanaka et al.l (2009) (T09 hereafter) subsequently pre- 
sented an exhaustive analysis of SN 2008D using a set of 
progenitor models and hydrodynamics/nucleosynthesis 
calculations. They found a very good agreement with 
the LC and the spectra for a model with Mej = 5.3 ± 1 
Mq, Ek=Q± 2.5 foe, and M{^^m) = 0.07 Mq. 

Early UV/optical emission is expected to occur after 
the arrival of the shock wave at the surface of the progen- 
itor (shock break-out) and before the re-brightening due 
to the decay of radioactive material. This emission is a 
consequence of the nearly adiabatic cooling due to the ex- 
pansion of the outermost layers of the ejecta. The obser- 
vations of the shock break-out and the following emission 
provide very valuable information about the structure of 
the star previous to the explosion. The duration of this 
early phase depends strongly on the size of the progeni- 
tor. For compact Wolf-Rayet stars, which are proposed 
progenitors of SNe lb and Ic, a duration of a few days is 
expected. Consequently, catching a SN during this phase 
is quite a challenge and only for a handful of objects has 
this been possible. 

Naturally, the origin of the earliest part of the 
LC of SN 2008D has been associated in the litera- 
ture to the adiabatic cooling of the outer layers of 
the eiecta (ISodcrb crg et al. 2008; Modjaz et al. 2009; 
iRabinak fc Waxmanl 120111) . However, different au- 
thors arrived at different conclusions regarding the pro- 
:enitor radius. Som e suggested a value of « 1 Rq 
Soderberg et al.ll2008HRabinak fc W axman"2011') while 
others proposed a larger value of « 9i?0 (Modia z et al.l 
12009V Che valier fc Franssonl l2008') . In ah these cases the 
estimations were bas ed on analytic models for the early 
emission. Note that iMazzah etall ()2008[ ) and T09 per- 
formed numerical simulations, but they did not attempt 
a hydrodynamical modeling of the early emission so as 
to help discriminate between both possibilities. 

Given the unique information about the progenitor 
structure that is provided by the early emission we 



have performed new hydrodynamical calculations for 
SN 2008D, focusing on this phase. Surprisingly, we found 
that our hydrodynamical models are not consistent with 
the cooling-phase explanation given in the literature for 
this object. Hydrodynamical models predict a much 
larger contrast between the luminosity minimum which 
occurs after the breakout, and the luminosity peak due to 
radioactive decay than what was ob served for SN 2008D . 
This has also been noted recently bv lDessart et al.l ()2011[ ) 
using a set of hydrodynamical and spectral calculations 
of SNe lb and Ic. 

Here we present the first radiation-hydrodynamical 
models for the early emission of SN 2008D and propose 
an alternative explanation for this emission based on a 
double-peaked ^^Ni distribution. The data and hydrody- 
namic code employed are described in § [2j A comparison 
with previous hydrodynamical results is shown in § 13.11 
Our proposed model for SN 2008D is presented in § 13.21 
In § 0] we compare our model with the analytic models 
used in the literature, and discuss alternative explana- 
tions for the early emission. Finally, our conclusions are 
presented in § [5l 

2. DATA AND MODEL 
2.1. Observational Material 

The bo lometric LC (Lbol) of SN 2008D was calculated 
by Modi az et aP ([2009) using UVWlBVRr'Ii'JHK, 
broad-band photometry for i < 31 days, referred to the 
onset of the XRT, i.e., JD = 2454475.06 (in what fol- 
lows, all times will be referred to this moment). Here we 
adopt the estimates of Lbol from black-b ody fits to the 
broadband photometry, which according to lModiaz et al.l 
(|200i l) are more accurate at early times than direct in- 
tegration of the observed flux. At later times we include 
bolometric calculations provided by T09 (see their Ap- 
pendix 1) using optical and near-infrared (NIR) data ob- 
tained with the MAGNUM telescope and the Himalayan 
Chandra telescope. Also available in the literature are 
two earlier data points, one at 0.14 days observed with 
the Swift telescope in the UVW2, UVWl and U fihers 
(jSoderberg et al.ll2008D. and anoth er at 0.44 days in the 
BVRI bands (jMazzah et al.ll2008D . The integrated flux 
for each of these observations was calculated respectively 
by the authors. Here we adopt, for t — 0.44 days, the 
sum of the luminosity in the BVRI bands, plus the con- 
tribution of the UV as estimated from the earlier Swift 
observations. The uncertainty of this point was assumed 
to be the sum in quadrature of the uncertainties of both 
contributions. We denote this earliest point with a dif- 
ferent symbol (square) in all figures where it is included 
in order to indicate its different origin. The distance and 
total reddening assumed in the calculations ar e d = 32 
Mpc, and E{B -V) = 0.65 mag (MazzaU et aLll2008[ ). 

The bolometric LC reveals a double-peaked shape 
that is also observed in the broad -band photometry 
(jSoderberg et al.|[200l [Mazzah et al.r 2008). We assume 
throughout this paper that the earliest data point is con- 
fident and therefore the early peak shape of the LC is 
robust. The first peak occurred at about 1 day with 
a luminosity of L = 1. x 10''^ erg s~^, and the sec- 
ond or main maximum happened at sal 9 days with a 
luminosity of L = 1.62 x 10^^ erg s~^, corresponding 
to Mboi = -16.8 mag. Therefore, SN 2008D had a nor- 
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mal peak luminosity, and a rise time at the long end of 
the observed range for SNe lb and Ic (|Richardson et alj 
120061 iDrout et al1l2012D . The contrast of luminosity be- 
tween the main peak and the dip that occurs before the 
re-brightening of the LC due to radioactive material is 
only of 0.35 dex. 

Apart from the bolometric LC, photospheric veloci- 
ties and color temperatures are also compared with our 
models. We employ He-line and photospheric velocities 
derived from the spectral modeling by T09, and color 
temperatures calculated by Modiaz ct al. (2009). 

2.2. Light Curve Models 

Synthetic LCs were calculated using the spherical, 
Lagrangian, LTE, h ydrodynamical code described by 
iBersten et all (|2011[ ) (BBHll, hereafter). The code 
solves radiation transport in the flux-limited diffusion 
approximation, including 7-ray transfer in gray approx- 
imation. Any distribution of the radioactive material is 
allowed and the energy deposition due to radioactive de- 
cay is computed in each layer of the entire ejecta. The 
ionization structure is determined by solving the Saha 
equation taking into account the most relevant elements 
in the progenitor structure. The Ross eland mean opac- 
ity is calculated using OPAL tables (jlglesias fc Roger j 
119961 and references therein), and the empirical relation 
used in T09, which is based on electro n scattering opacity 
as derived from t he TOPS database (jMagee et al.lll995l : 
iDeng et aLll2005[ ). and including effects of lines. The ex- 
plosion is simulated by injecting near the center of the 
object a certain amount of energy, in a thermal form, 
during a short interval. The code does not explicitly 
solve for the explosive nucleosynthesis produced during 
the shock propagation, but it implicitly takes this into 
consideration in the chemical composition assumed for 
our initial models. 

As pre-supernovae models, stellar evolu ti on ca lcula- 
tions of He stars by iNomoto fc Hashimotol (|1988f ) were 
adopted. T09 studied five different He star models with 
masses of 4, 6, 8, 10, and 16 Mq, using a variety of 
explosion parameters (see their Table 1). They found 
that the models with He mass of 6 and 8 Mq (He6 and 
He8, respectively) were most consistent with the spec- 
tra and LCs of SN 2008D. They derived a kinetic en- 
ergy of Ek = (6.0 ± 2.5) X 10^^ erg, an ejecta mass of 
= 5.3±1.0 Mq, a mass ofM{^^m) = 0.07 Mq, 
and a compact-remnant mass in the range of 1.6 — 1.8 

Mq. 

We have based the current analysis on the results of 
T09, and so we calculated the hydrodynamics for these 
two optimal models (He6 and He8), which correspond to 
main-sequence masses of 2 and 25 Mq, as deriv e d from 
the Mms - Ma relation of lSugimoto fc Nomotd (|1980( ). 
The chemical abundance distribution left by the explo- 
sive nucleosynthesis is assumed as a pre-explosion con- 
dition (see T09 for more details). One characteristic of 
these nucleosynthesis calculations in spherical symme- 
try is that '^^Ni is confined in the innermost layers of 
the ejecta, which makes it very difficult to explain the 
timescale of the rising part of the LC normally observed 
for SNe lb and Ic. Instead, a more extensive mixing of 
^^Ni, associated to multi-dimensional effects, was used 
by T09, and we also adopted it here (see ? 13.11 for more 



details). Note that mixing of ^^Ni to large radius of 
the ejecta of CCSNe has been suc cessfully produced in 
recent 3D numerica l calculations (iHaminer et al.' '201d; 
iJoggerst et al.l 120101 ). In addition, Pessart et al. (2013) 
have claimed that all SNe classified as lb require efficient 
mixing of ^^Ni in the helium-rich layers. 

3. HYDRODYNAMIC MODEL OF SN 2008D 

We present here our hydrodynamic calculations for 
SN 2008D using the He6 and He8 initial models. In 
§ 13.11 we compare our results with those of T09 by fo- 
cusing on the second peak of the bolometric LC (t > 5 
days), which determines the global pro perties of the SN. 
The early emission is analyzed in § 13. 2[ where we present 
a model that consistently reproduces the the first and 
second peaks of the LC. 

3.1. Comparison with T09 

As it is well known, the shape of the LC depends on 
the explosion kinetic energy Ek, the ejected mass Moj, 
and the mass and distribution of ^^Ni. A nalytic expres- 
sions for these dependencies were given bv lArnettl (|l982'). 

where the width of the LC peak is tlc oc M^j'^ Ej^^^'^ , 

and the peak luminosity is Lpeak oc AfNiTj^Q. There- 
fore, it is possible to estimate explosion parameters by 
comparing models and observed bolometric LCs. How- 
ever, various combinations of Moj and Ek can fit the 
LC, and spectra modeling is needed in order to break 
the degeneracy of the parameters. This type of analysis 
was done by T09 and they found a very good agreement 
between models and observations using the He6 and He8 
initial models with the following explosion parameters: 
for He6, Ek = 3.7 foe, a cut mass of Mcut = 1.6 Mq, 
and Mm = 0.065 Mq; and for He8, Ek = 8.4 foe, 
Mcut = 1.8 Mq, and Mm = 0.07 Mq. A constant ^"^Ni 
distribution up to 7000 (9000) km s~^ for model He6 
(He8) was assumed to account for the rise time to the 
main peak. 

We have used the code of BBHll to calculate bolomet- 
ric LCs and photospheric evolution for the same models 
(He6 and He8) and explosion parameters as those of T09. 
As opposed to T09, our calculations solve the hydrody- 
namics coupled to the radiative transfer, allowing us to 
model consistently the earliest phases of the SN evolu- 
tion. Figure [1] shows a comparison of our results with 
those of T09, along with the observed bolometric LC 
of SN 2008D. The photospheric velocities are compared 
in Figure O where we include the He lines velocities as 
well as the photospheric velocity estimated from spectra 
modeling by T09. From these figures we see that the 
agreement between both models is reasonable, consider- 
ing the differences in the calculation methods. The LCs 
look remarkably similar for times t > 5 days. It is not 
striking that the largest differences appear at the earliest 
epochs, since the code in BBHll self-consistently calcu- 
lates the shock wave propagation, the breakout and the 
later evolution, while in the calculations of T09, the hy- 
drodynamics and the radiative transport were computed 
with different codes, switching from one code to other 
when the homologous expansion was achieved. There is 
a small systematic difference between the photospheric 
velocities of both models. This is probably related to 
slight differences in the re-gridding of the initial model. 
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Fig. 1. — Comparison between the bolometric LC for models 
He6 (blue lines) and He8 (red lines) calculated in this work (solid 
lines) and by T09 (dashed lines). Both works adopted the same 
progenitor models and physical parameters. The di ffere nces are 
only related with the hydrodynamic code used (see § 13.11 for more 
details). The observed bolometric LC of SN 2008D is also included 
(cyan dots). 
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Fig. 2. — Evolution of the photospheric velocity for models He6 
(blue lines) and He8 (red lines) calculated in this paper (solid lines) 
and by T09 (dashed lines). For comparison, observed He line ve- 
locities (circle dots) and photospheric velocities estimated through 
spectra modeling by T09 are also shown. 



which produces a smah difference in the actual total mass 
of the progenitor (of w 0.4 Mq). beginfigure 

Both our calculations and those of T09 provide very 
good matches to the observations around the main peak. 
Note that model He8 is the one that best represents the 
bolometric LC, while model He6 gives a better fit to the 
velocities, specially the photospheric velocities. Some in- 
termediate model between these two seems to be the 
most plausible progenitor, as suggested in T09. How- 



ever, at times before 8 days, the models predict much 
lower luminosity than the observed one. As described in 
§ 13. 2[ this can be overcome by including a small amount 
of ^^Ni in the outermost layers of the ejecta. 

3.2. A Htgh- Velocity blob of ""^Ni 

It is clear from Figure [1] that the models presented in 
ij 13. II cannot explain the first peak shown by the observa- 
tions. The difference in luminosity between the hydrody- 
namical models and the observations at early times are 
larger than 0.5 dex, which is much greater than any possi- 
ble uncertainty related with the calculations and observa- 
tions. The models predict a brightness contrast between 
the dip and the main peak that is larger than 0.9 dex, 
while the observations show a contrast of only 0.35 dex. 
Although the early behavior of SN 2008D was attributed 
in the literature to the cooling of the outer envelope of the 
ejecta after being heated by the shock wave, our calcu- 
lation suggests that this cooling happened very quickly, 
and in less than 1 day most of the energy deposited by 
the shock was degraded. Therefore, unless we assume a 
different structure for the progenitor — i.e., lager radius 
and/or more massive external envelope — , some source of 
extra energy in the outer layers of the ejecta is required 
in order to reproduce the first peak of the LC. Here we 
study this possibility and leave the discussion of other 
alternatives for § 21 

We artificially placed some ^^Ni in the outer layers of 
the ejecta as the source of extra energy needed to ex- 
plain the first peak. This material may have been car- 
ried by a jet-like phenomenon produced during an as- 
ph erical explosion. The presence of a jet was proposed 
by IMazzali et al.l ()2008D to explain the XRT associated 
with SN 2008D, as well as the broad-line appearance of 
the spectrum at t < 3 — 5 days. T he double-pe aked oxy- 
gen lines detected in late spectra (iMalesani et al. 2009| ; 
iModiaz et all 120091 ; iTanaka et all l2009f ) provide addi- 
tional evidence of the asphericity during th e explosion. A 
spectr opolarimetry study of SN 2008D by iMaund et al.l 
(j2009| ) suggest that a jet was produced but that it stalled 
in the C+0 core. Maund et al. base this conclusion 
mainly on the observed low degree of O I line polariza- 
tion. However, for this SN the O I A 7774 line is so weak 
that one cannot expect any strong polarization associ- 
ated with it. Weak lines do n ot produce high pola riza- 
tion, as shown, for example, inlTan aka et al.l ()2012[ ). 

A double-peaked ^^Ni distribution was previously 
suggested t o model the double-peaked light curve o f 
SN 2005bf (iTominaga et al.l l2005t IfblateUi et al.l l2()0fil) . 
The presence of an unobserved jet in SN 2005bf was spec- 
ulated as being responsible for the anomalous ^^Ni dis- 
tribution. Despite the overall observational differences 
between SN 2005bf and SN 2008D, some similarities in 
the polarization pro perties of both objects were found 
(jMaund et al.l 120091 ) that could indicate similarities in 
the explosion geometry. 

In Figure |3| we show the bolometric LC for our model 
with ^^Ni in the outer layers (solid line) compared with 
the observations. This "optimal model" is similar to the 
He8 model presented in § 13.11 (dashed line in the figure) 
but with an extra amount of 0.01 Mq of ^^Ni in the out- 
ermost layers of the ejecta (at v > 20, 000 km s""'^) added 
to explain the early emission at t < 5 days. In addition. 
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Fig. 3. — Bolometric LC of our optimal model (solid line) com- 
pared with the observations of SN 2008D. This model can con- 
sistently reproduce the two observed peaks. For comparison, two 
models without external ^''Ni and slightly different internal ^''Ni 
distributions are shown. A model with ^^Ni mixed up to 9000 km 
s~^ as assumed in TOO is shown in dashed line, and a model with 
the same internal ^^Ni distribution used in the optimal model, i.e. 
up to Si 10000 km s~^ is shown in dotted line. Note that both 
models without the external ^®Ni show a difference with the ob- 
servations larger than 0.5 dex at early times. 



the optimal model has a slightly extended distribution of 
internal ^^Ni up to w 10000 km s~^ (as compared with 
9000 km s""'^ that was used in the model of T09) to im- 
prove the agreement with the data during the rise to the 
main peak. For clarity, in Figure |3] we show two mod- 
els that lack external ^^Ni, with internal distributions up 
to 9000 and 10000 km s~^ (dashed and dotted lines, re- 
spectively). In the rest of this paper the model without 
external ^^Ni will be the latter one and will be referred 
to as He8. 

The agreement between the optimal model and the ob- 
servations is excellent. The assumption of external ^^Ni 
allowed us to reproduce the first and second peaks con- 
sistently. The luminosity contrast between dip and main 
peak is now similar to that shown by the observations. 

In order to obtain our optimal model, several ^^Ni dis- 
tributions were explored. In Figure |4] we schematically 
show a double-peaked ^^Ni distribution where the inter- 
nal and external components were assumed to be step- 
like functions within some internal and external veloci- 
ties. This distribution can be characterized with three 
parameters: the mass of external ^^Ni, M(^^Ni)out, the 
minimum velocity for the external ^^Ni, u(^^Ni)out, and 
the maximum velocity for the internal ^^Ni, u(^^Ni)in. 
The effect of the variation of these parameters on the 
early LC is shown in Figure [3 We see that (1) larger 
M(5*'Ni)out produces a more luminous first peak, (2) 
higher u(^^Ni)out produces an earlier first peak, and 
(3) lower w(^^Ni)in translates to deeper and later min- 
imum, or later rise to the second peak. Specifically, we 
have adopted the following values for our optimal model: 
M(56Ni)out = O.OIMq, w(56Ni)out = 22000 km s^^, and 
i,(56Ni)i„ = 10000 km s-i. 

Although we do not provide a detailed calculation to 
account for the occurrence of such ^^Ni distribution, it is 
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Fig. 4. — Schematic doubly-peaked ^^Ni distribution charac- 
terized by the following parameters: the mass of internal ^^Ni, 
Af(56Ni)in, the mass of external ^''Ni, M(5^Ni)out, the minimum 
velocity for the external ^^Ni, f(^®Ni)out and the maximum veloc- 
ity for the internal ^SNi, D(5<3Ni)in. 
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Fig. 5. — Sensitivity of the early LC on the distribution of ^®Ni. 
Three different values of Af(f '5Ni)out , ■«(^'^Ni)out and t;(^'*Ni)in 
have been used (see Figure |4| . The arrows indicate the change 
produced in the LC by the increase of each parameter. A model 
without ^^Ni is shown for comparison. Time is plotted in logarith- 
mic scale to show more clearly the behavior of the early LC. 



remarkable that it allows us to obtain a very good fit to 
the observation and also to explore its effect on the early 
LC. We also explored the possibility of using smoother 
functions for the ^^Ni distribution, but we found that 
they failed to reproduce the two peaks shown in the ob- 
servations, and instead they tended to produce a plateau- 
like shape. Given the large uncertainties in the early ob- 
servations, we emphasize that the critical feature of our 
model is to have some '^^Ni in the outer layers but not 
the exact amount and shape of its distribution. 

Another important parameter to compare with the ob- 
servations is the temperature. Figure IH] shows the effec- 
tive (solid lines) and color (dashed lines) temperatures 
for the optimal model (red) and the model without exter- 
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Fig. 6. — Color (dashed lines) and effective temperature (solid 
lines) evolution for our optimal model (red) and for the same model 
without external ^^Ni (blue). The color temperature of SN 2008D 
calculated by Modiaz et al.! I|2009l ) using broad-band photometry 
is shown for comparison (dots). 



^^Ni. Therefore, we do not present a comparison between 
the velocities of the optimal model and observations be- 
cause this is essentially the same that we presented in 
§[3l]and Figure El 

One could expect that the presence of ^^Ni and asso- 
ciated iron-group elements in the outer layers can intro- 
duce lines and increase the blanketing, thus affecting the 
observed spectrum. However, if the iron-rich material is 
confined to a small solid angle, as in the case of a jet, 
the effect will be diluted, as explained bv iDessart et alJ 
iOlJ). 

As our simulations are carried out in one dimension 
we do not have the ability of reproducing the actual 
structure of a jet-like distribution. Nevertheless, in the 
proposed scenario only a small fraction of the mass is in- 
volved in the jet itself, thus not affecting the global explo- 
sion dyna mics. The low level of continuum polarization 
found by fiMaund et al.l 12009) indicates that departures 
from spherical symmetry should be small. This allows 
us to treat the problem approximately in spherical sym- 
metry. Calculations performed in higher dimensions are 
required to produce a self-consistent model that would 
allow to test the proposed jet and aspherical ^^Ni distri- 
bution. 



nal ^^Ni (blue), compar ed with the observe d color tem- 
perature of SN 2008D (jModiaz et all [20091) . The color 
temperature (Tc) gives information about the contin- 
uum spectral energy distribution. To calculate Tq from 
the models it is needed to know the layer in the ejecta 
where the spectrum is formed. An approximate estima- 
tion of Tc is gi ven by the temperature at the "thermal- 
ization" depth (|Ensman fc Burrowslll993) . below which 
the gas and radiation field are in equilibrium. Following 
lEnsman fc Burrowsl (|1992| ). we calculated the "thermal- 
ization" depth as the layer where SrabsTsct ~ 1- Here 
Tsct is the optical depth for scattering and Tabs the opti- 
cal depth for absorption determined using /tabs — k, — ^sct 
where k is the Rosseland mean opacity calculated with 
OPAL tables, and Ksct is the scattering opacity calcu- 
lated by solving the Saha equations. The effective tem- 
perature is defined as T^^ = L/{4tt aR^^), where i?ph is 
the radial coordinate at the photosphere position defined 
as the layer where t — 1 and L is the radiative luminos- 
ity plus the luminosity of ^^Ni decay deposited above the 
pliotosphereQ- From Figure |6] we see, as expected, that 
the model without external ^^Ni has higher color than ef- 
fective temperature, but neither of them provides a satis- 
factory agreement with the data. On the other hand, the 
temperatures of the optimal model compare acceptably 
well with the data, but unlike what is expected, the ef- 
fective temperature is higher than the color temperature 
for t ^ 1 day. This can be understood from our def- 
inition of the effective temperature, which includes the 
extra luminosity of ^^Ni decay above the photosphere. 

Unlike luminosity and temperature, the photospheric 
velocity is almost unaffected by the existence of external 

^ As stated in Section l2.2l we have solved the gamma-ray trans- 
fer which allow us to compute the actual deposition of energy of 
gamma-ray in the entire ejecta (in particular, above the photo- 
sphere). Such deposition can be appreciably lower than the total 
available radioactive energy release, specially when the envelope 
becomes very dilute. 



4. ALTERNATIVE APPROACHES 
4.1. Cooling of the Shocked Envelope 

In the literature, the early emission of the LC of 
SN 2008D was explained as a consequence of the cool- 
ing of the outer stellar envelope following the passage 
the shock through the st ar and its subsequent break- 
out. Analytic models bv |^yaxman etall (|2007D (W07) 
and 'Chevalier fc Franssorf ('2008 D (CF08) were used in 
Sodcrbcrg ct al. (2008) and Modi az et all (I2009D to com- 
pare with the early data of SN 2008D. These models 
describe the emission of the outer layer of the ejecta 
assuming: (1) pre-explosion density p oc (1 — r/i?)", 
which is valid while the photosphere is in the outer shock- 
accelerated part of the ejecta as long as the mass above 
of the photosphere is less than 0.1 Mq, (2) self-similar 
solution once the supernova reached the state of free ex- 
pansion to determine the post breakout density and ve- 
locity, and (3) constant opacity. The luminosity was cal- 
culated in a different way in each model. While W07 did 
not take into account the radiative diffusion assuming a 
strictly adiabatic expansion, CF08 considered the motion 
of a diffusion wave through the ejecta. Therefore these 
models are valid after free expansion is achieved and un- 
til more or less the onset of recombination, when the 
photosphere begins t o rece de into the ejecta. Recently, 
Rabinak fc Waxma3 (|2011[ ) (RWll) improved the model 
of WO 7 by including the effect of recombination on the 
opacity. They also corrected a typographical error found 
in a coefficient of the equation (19) of W07. 

The analytic models mentioned above provide expres- 
sions for the time evolution of the luminosity, photo- 
spheric radius and temperature as a function of ejecta 
mass (Mcj), ejecta kinetic energy (Ek) and progenitor 
radius (i?), besides other parameters that depend on the 
structure of the progenitor. However, the dependence 
of luminosity and temperature on R is stronger than on 
the other parameters. It is linear for the luminosity and 
cx i?^/^ for the temperature. 
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Fig. 7. — Comparison between the bolometric luminosity of the 
analytic and hydrodynamical models. The data are shown with 
cyan dots. Three different hydrodynamical calculations are shown: 
the optimal model (solid line), the same model but without exter- 
nal ^"Ni (dashed line), and a model assuming a constant opacity 
equal to the electron-scattering opacity for pure helium material, 
Kg- = 0.2 g cm~^, (dotted line). The analytic models of CF08 
(black) and RWll (blue) for the same explosion parameters as in 
our hydrodynamic simulations are also shown (solid line). With 
dashed lines we show the same analytic models for an alternative 
set of physical parameters that include a larger initial radius. Note 
that a larger radius provides a better fit to the observations, with 
the exception of the earlies t point . The analytic prediction of the 
binary interaction model of lKasenl 12013) is also included (see ij |4.3l 
for additional information). 



In Figure [7] we show a comparison between our hydro- 
dynamical calculations and the analytic predictions of 
CF08 and RWll. Three different hydrodynamical calcu- 
lations are shown: the optimal model, the same model 
without external ^^Ni, and a model with constant opacity 
equal to the electron-scattering opacity for pure helium, 
Kg- = 0.2 g cm^^. For the analytic models, the corrected 
expression for constant opacity from RWll [their equa- 
tion (13)] was used here instead of equation (19) of W07. 
Two sets of explosion parameters for models CF08 and 
RWll are shown. First, we focus on the behavior of the 
analytic prediction with the same explosion parameters 
as in our hydrodynamic simulations (i.e., Ek — 8.4 foe, 
Mej = 6.2 Mq, and R = IARq). The analytic models for 
compact stellar structures assume an initial density pro- 
file with index n = 3, which is similar to the shape of the 
initial structure from stellar evolution that we adopted 
for model He8, as shown in the inset of Figure HI As the 
proportionality constant between p and r is not given in 
the analytic models, we have scaled the density of these 
models to match the value of the He8 model at 0.1 Mq 
inside the star, i.e. where the analytic models are valid. 

From Figure [7] it is clear that (1) the analytics mod- 
els predict a higher luminosity than the hydrodynamic 
model without external '^^Ni, but they are still not con- 
sistent with the observations, and (2) the slope of the 
analytic and hydro models are quite similar until t w 0.5 
d, when helium starts to recombine and the luminos- 
ity of the hydro model enters a nearly plateau phase for 
approximately five days. After that time, the heating 



by ^^Ni produces a re-brightening of the LC. Note that 
the presence of a post-breakout plateau for Wolf-Rayet 
progenitors of different radii an d masses has recently 
been reported in simulations by iDess art et al.l (pOll). 
In agreement with our simulations, they found that the 
typical differences between the luminosity of the post- 
breakout plateau (1 — 5 x IQ^Lq) and the main peak 
are much larger than the observed value of « 0.3 dex 
for SN 2008D. Our hydrodynamical models show that 
the hypothesis of constant opacity breaks at i « 0.5 d, 
when He begins to recombine and the photosphere re- 
cedes into the ejecta. By t « 1.5 d, the mass above 
of the photosphere is larger than 0.1 Mq and the as- 
sumptions of the analytical models are no longer valid. 
Therefore, our simulations establish a limit of about 1.5 
days for the validity of the analytic expressions. This is 
in close concordance with the range of validity suggested 
by equations (16) and (17) of RWll, which for the val- 
ues of Ek, Mej and R used here, give a range of up to 2 
days. Note that for this range of time there are only two 
data points available, including the earliest point which 
seems not to follow the analytic predictions. This makes 
the conclusions derived from the analytic models more 
dubious. 

Finally, as we noted before, there is a difference in lumi- 
nosity between our models and the analytic models even 
for times earlier than 0.5 day. The differences may be re- 
lated with the density structure of the outermost layers of 
the ejecta used in each calculation. A direct comparison 
of initial density structures is not possible because the 
scale of the relation between p and r is not provided in 
the analytic calculations. We can however make a quan- 
titative comparison of post-shock breakout density struc- 
tures. For the analytic models, this has the form p oc w". 
Equation (1) of CF08 evaluated at E ^ 8.4 foe, M = 6.4 
Mq, and t = 1.1 days gives p = 7.9 x iQ^y-'^^o.is g ^m-s^ 
where v is expressed in 10® cm s^^. For our density pro- 
file at i = 1 . 1 days we find a similar exponent but a very 
different proportionality constant: p = 1.9 x 10^ 
g cm"-^. This difference can be the reason for the dis- 
crepancy in luminosity between analytic and numerical 
models. The difference in density may occur because 
the calibration of the post-explosion density profile used 
in CF08 was based on previous hydrodynamical calcu- 
lations applied to a blue supergiant structure, useful to 
model SN 1987A, instead of a Wolf-Rayet progenitor as 
required here. 

4.2. Different Progenitor Structure 

A larger value of the progenitor radius of 9 Rq was sug- 
gested by CF08 in order to explain the XRT of SN 2008D 
as completely therm al emission. On the other hand, 
iModjaz et al.l ()2009f ) found a similar radius by fitting the 
CF08 relations to the early data, and a ssuming Ek = 2 
foe an d Mej = 5 Mq, as suggested bv iSoderberg et aLl 
((200l) . In Figure El we include the analytic models for 
these alternative physical parameters. In this case, the 
analytic models reproduce better the early luminosity 
with the exception of the earliest data point. 

The improvement in the fit of the analytic models with 
larger radius suggests that adopting progenitors with 
large radii in our hydrodynamical simulations may al- 
low us to reproduce the early data without resourcing to 
external ^^Ni. We thus attached several envelopes in hy- 
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Fig. 8. — Initial density distributions as a function of radius for 
model He8 {R = 1.4 /Jq; red line) and for models with different 
initial radii. Thick lines represent models whose variations in ra- 
dius are accomplished by attaching essentially mass-less (< 0.01 
Mq) envelopes to the He-rich layer of HeS model while thin lines 
show models with 9 Rq and massive envelopes attached at differ- 
ent points of mass inside the HeS model, as indicated in the upper 
right part of the figure. Inset: blow out of the outermost layers 
for model HeS (R = 1.4 Rq; red line) compared with the shape 
assumed for the analytic models: p oc (1 — r/R)", with n = 3 and 
R = 1.4 /Jq (dashed black line). The location of the limit of va- 
lidity of the analytic models (0.1 Mq inside the star) is indicated 
with vertical lines. 




t [days] 



Fig. 9. — Bolometric LC for models with the same explosion 
energy as the optimal model, but different initial radii. Models 
whose variations in radius are accomplished by attaching essen- 
tially mass-less (< 0.01 Mq) envelopes to the He-rich layer of HeS 
model are shown with thick lines. Note that in this case, larger 
radii produce higher early luminosity for t < 6 days, although not 
high enough to explain the first peak of SN 200SD. Thin lines rep- 
resent models with 9 Rq but with massive envelopes attached at 
different points of mass inside the HeS model, as indicated in the 
labels of the figure. These models, especially Menv = 5 and 6 Mq 
give a reasonable fit to the early observations with the exception 
of the earliest data point. 



drostatic and thermodynamical equilibrium to the He- 
rich layer of our He8 model. The attached envelopes 
were integrated inward assuming a fixed stellar radius 
and varying the effective temperature and envelope mass 
so that the condition of continuity at the fitting point 
was achieved for mass, luminosity, pressure, and temper- 
ature. We also verified that the bottom of the envelope 
is cool enough to prevent nuclear reactions from devel- 
oping in that region. This allowed us to generate a set 
of progenitors with radii of 9, 20, 50 and 100 i?©. 

The resulting bolometric LCs are shown with thick 
lines in Figure IHl where the other physical parameters 
are the same as previously adopted. From the figure it is 
clear that at t < 5 days models with larger radii produce 
slower cooling of the outer layers and higher luminosity, 
though not as large as the one predicted by the analyti- 
cal models. After that time, all models look remarkably 
similar. Note that models start to converge at t w 1.5 
days when the analytic expressions become invalid. Fi- 
nally we see that even with a radius as large as 100 Rq, 
we could not satisfactorily reproduce the early LC. 

The analysis above was performed using nearly mass- 
less envelopes (< 0.01 Mq). The initial density profiles 
as a function of radius for these models are shown with 
thick lines in Figure [H The shape of the density in the 
outermost layers not to follow the one used for the ana- 
lytical models. This is reflected in a different, smoother, 
dependence of the luminosity on the progenitor radius 
than the prediction of the analytic models. 

Alternatively, we calculated initial models with 9 Rq 
but with a substantial modification of the initial density 
distribution with respect to the He8 model. Specifically, 



we removed 4, 5 and 6 Mq of the original He8 model 
and attached new massive envelopes preserving the total 
mass (8 Mq), out to a radius of 9 Rq. These structures 
are shown in Figure ITOl as a function of mass. The mass- 
less envelope model for 9 Rq is also included for compar- 
ison. Figure [S] shows the LCs resulting from these models 
(thin lines). The early observations can be reproduced 
reasonably well with these type of structures, excluding 
the earliest data point. In these cases, the shape of the 
density profile as a function of radius in the outermost 
layers (thin lines in Figure [8|) are closer to the function 
assumed for the analytic models. 

This means that, without considering the earliest data 
point, models with the larger radius and a modified den- 
sity structure than those predicted by stellar evolution 
calculations can give a reasonable explanation of the 
early LC of SN 2008D as the cooling expansion of the 
outer envelope. However, the massive-envelope models 
give a much poorer fit to the LC around the main peak, 
and they pose the additional problem of finding a phys- 
ical explanation to justify a density profile that is dif- 
ferent from those predicted by stellar evolution models. 
Note that some variations in the assumed pre-explosion 
density profile can be caused by effects of rotation, which 
were not included in our initial models. Nevertheless, the 
expected differences are much smaller than those needed 
to fit the early data of SN 2008D. 

It is important remark that r ecent binary stellar evo- 
lution calculation presented by lYoon et al.l (|2010D pre- 
dicted a larger progenitor radius than that of 1.4 Rq 
obtained for our HeS model based on single stellar evo- 
lution. However, the radius comprising 95% of the mass 
is less than 1 Rq , with the exception of the models that 
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Fig. 10. — Initial density distributions as a function of mass for 
models with 9 Rq formed from model He8 but with massive en- 
velopes attached at different points of mass, as indicated in the 
labels. 



contain some hydrogen for whi ch this radius can be as 
large as 5 Rq (see Table 2 of lYoon et aT1l20Tol ). For 
SN 2008D, the presence of a thin H envelope was firmly 
ruled out from spectroscopic analysis, e.g. T09 estimated 
a very low upper limit for the hydrogen mass fraction of 
5 X 10"'* Mq. Therefore, the binary models are similar 
to our low-mass envelope models and are not consistent 
with the early observations of SN 2008D. This also be- 
comes clear from the results present ed bvlDessart et all 
(|2011h based on the binary models oflYoon et al.l (I2010D. 
All the LCs presented by iDessart et al.l ()2011f ) have a 
post-breakout luminosity that is much smaller than the 
one observed for SN 2008D, as noted by the authors and 
according also to our calculations. 



4.3. Binary interaction model 

As we noted previously, IDessart et al.l ()2011D also 
found that the observed luminosity of SN 2008D at early 
times was much larger than the one estimated by their 
models. Their calculations were done with a different 
hydrodynamical code and evolutionary initial models, as 
compared with our study. They considered unlikely that 
light contamination from the host galaxy, or any light 
scattered by the CSM or pre-SN mass loss could produce 
the large post-breakout luminosity. Moreover, CF08 es- 
timated that the mass-loss rate of the progenitor was 
too low for th e CSM to be optically thick. Alternatively, 
IDessart et al.l (j2011l ) suggested that the large observed 
luminosity could be due to the collision of the SN ejecta 
with a companion star in a binary system. However, they 
did not perform any detailed test of this scenario. Mod- 
els that explain the early enhancement of the luminosit y 
due to binary colhsion were proposed bv iKaseiil ()2010[ ). 
These models depend on several parameters, such as bi- 
nary separation, mass of the ejecta, shock velocity and 
viewing angle. An analytic expression for the luminos- 
ity as a function of time for a viewin g angle of 45° is 
given in Equation (22) of iKasenI (|20100 . with a stronger 
dependence on the shock velocity and the binary separa- 
tion (a) than on the ejected mass and the electron opac- 



ity. The shock velocity for SN 2008D can be estimated 
from the ejecta mass and explosion energy, assuming that 
Kg- = 0.2 g cm~'^ for fully ionized electron-scattering of 
pure helium, and thus leaving a as the only free param- 
eter of the problem. Figure [7] shows an example of the 
binary collision model for a shock velocity derived using 
an explosion energy of 8.4 foe and an ejecta mass of 6.2 
Mq as adopted in this paper. For this set of parameters 
the orbital separation required to reproduce the early 
luminosity is a = 1.8 i?Q. This value can be slightly in- 
creased to a = 2.8 Rq if we assume values of E = 3.7 foe 
and Af cj = 4.4 Mq appropriate for the Hc6 model (see 
section |3T|) ■ This very small constraint found for the bi- 
nary separation may pose a problem for the binary colli- 
sion scenario. Note that the typical binary separation of 
Galactic Wolf-Rayet stars is > IORq, with very few ex- 
ceptions (jvan der H ucht 2001). The situation might be 
remedie d by assumin g smaller viewing angles (see Fig- 
ure 2 of IKasenI [2OT0h . but since we have no way to test 
this scenario, we will consider it a less likely possibility. 

5. CONCLUSIONS 

The early post-breakout emission of SN 2008D can- 
not be explained as the cooling of the outer layers of a 
normal Wolf-Rayet star after the passage of the shock 
wave, as suggested in the literature. Alternatively we 
have proposed a model which assumes a double-peaked 
^^Ni distribution to explain this early emission as an at- 
tractive solution to this problem. Specifically, an amount 
of 0.01 Mq of ^^Ni was located in the outermost layers of 
the ejecta with v > 20, 000 km — although the exact 
amount and distribution is subject to some uncertainties 
inherent to the data and the model itself. The assump- 
tion of external ^^Ni allowed us to reproduce very well 
both the early and late observations of SN 2008D. The 
presence of this high-velocity radioactive material may 
be caused by the formation of jets during the explosion. 
A multidimensional model would be required to confirm 
our suggestion in a self-consistent fashion. 

Hydrodynamical models applied to Wolf-Rayet star 
predict a much larger contrast (^0.9 dex) between the 
initial dip and the peak of the LC due to heating by 
than what is observed for SN 2008D (« 0.3 dex). 
Only with a substantial modification of the initial den- 
sity structure that is predicted by stellar evolution cal- 
culations and for a larger progenitor radius, 9 Kq, can 
the early emission be compatible with the cooling of the 
outer envelope. Even in this case the earliest observed 
point is not reproduced by the shock-cooling model and 
a much worse fit to the LC around the main peak is ob- 
tained. In addition, some physical explanation for such 
substantial changes in the initial density is needed to 
consider this as a possibility. 

We also analyzed the possibility that the early emis- 
sion was due to interaction of the ejecta with a binary 
companion. This was base d on the analy tic predictions 
for the luminosity given bv IKasenI ()2010l ). However, we 
found that the binary separation required to explain the 
early luminosity of SN 2008D is very small (< 3Rq), 
which poses a serious conflict for this interpretation. 

Furthermore, we note that only the double-peaked ^^Ni 
model can reproduce the earliest observed data point as 
well as the rest of the LC. Both the binary interaction 
model and the modified density structure predict a larger 
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luminosity at the earliest observed epoch. 

Comparing our hydrodynamical models with the ana- 
lytic models commonly adopted in the literature to re- 
produce the early emission (e.g.; CF08 and RWll) we 
found: (1) the range of applicability of the analytic for- 
mula is very restricted. An upper limit of w 1.5 days 
after the explosion was found. In this range of time only 
two data points of SN 2008D are available. (2) The an- 
alytic models were calibrated with structures that may 
not apply to Wolf-Rayet stars. This could be the reason 
for the very different value found for the proportional- 
ity constant of the relationship p oc for the post- 
explosion density profile as compared with the one used 
in CF08 (Equation 1). (3) A modification of the radius 
of the star, as commonly applied for the analytic mod- 
els, implies a different inner boundary condition for the 
density which cannot represent any stellar evolutionary 
solution. Therefore, conclusions based on these simple 
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models should be taken with caution. 

The analysis presented here shows the relevance of the 
early emission in our understanding of the progenitor 
structure previous to the explosion. We expect that cur- 
rent SN searches of increasing cadence will frequently de- 
tect more CCSNe in the early stages of their evolution. 
Then it will be possible to test whether the early behav- 
ior of SN 2008D is a peculiar to it or if it is a common 
feature of some type of SNe. 
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